Translation invariance of weak KAM solutions of 
the Newtonian n-body problem 

Ezequiel Maderna 
January 23, 2013 

Abstract 

We consider the Hamilton- Jacobi equation H{x, d^u) = c, where c > 0, 
of the classical A'^-body problem in some Euclidian space E of dimension 
at least two. The fixed points of the Lax-Oleinik semigroup are global 
viscosity solutions for the critical value of the constant (c = 0) , also called 
weak KAM solutions. We show that all these solutions are invariant under 
the action by translations of E in the space of configurations. We also show 
the existence of non-invariant solutions for the super-critical equations 
(OO). 

1 Introduction and results 
1.1 Preliminaries 

Let E be an Euclidean vector space of dimension fc > 2 in which N punctual 
masses mi , . . . , toat > are mutually attracted by the Newtonian gravitational 
forces. At each time, the position vectors of the N bodies determine a config- 
uration of the system x — (ri, . . . ,rjv) G E'^ . The Newtonian potential is the 
function U : E^ — >• (0, -l-oo] defined by 



U{x) = 

l<i<j<N 



niimj 



It is clear that U{x) < +oo if and only if a; is a configuration without collisions, 
meaning that there are no values of i < j such that = rj. Moreover, if 
we denote by O C E^ the open and dense subset of configurations without 
collisions, then the Newtonian potential is an analytical function on fi, and the 
law of gravitation is equivalent to the second order differential equation x = VU 
on f2, where the gradient is taken with respect to the mass inner product. 

As well as the equations of motion, the Hamiltonian and the Lagrangian of 
the system admit a synthetic expression in terms of the mass inner product. We 
recall that given two configurations x = (ri, . . . , r^r) and y — (si, . . . , sn), their 
mass inner product is given by 

N 

X ■ y ^ y^ rui < n, Si > . 
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With the natural identification TE^ ~ x E'^ , the Lagrangian function of 
the system reads 

Lix,v) = ^ \ vf + U{x) 

where | ?; | = (v ■ ?;)^/^ denotes the induced norm by the mass inner product on 
. Consequently, if p = (pi, . . . ,pAr) is in the dual space {E'^)* ~ (E*)^ , and 
I p I denotes the operator norm of p for the mass inner product, we have 

where \\pi\\ is the operator norm of pi € E* . Since the Hamiltonian H of the 
system is defined on each fiber of the cotangent bundle T*E^ as the convex 
dual of the Lagrangian, we obtain the expression 

H{x,p) = \\p\^ -U{x). 

Recall that here the Legendre transform u) = dvL{x,v) 6 T*E^ does not 
depends on the configuration x, hence it can be defined as a map C : E^ — >■ 
{E^y. Thus p = C{v) will mean that p{'w) = w • w for all w e E^ . It is easy to 
check the Young-Fenchel inequality, which assures that for all p £ (E^)* and 
V G E^ we have p{v) < H{x,p) + L{x, v), and that equality holds if and only if 
P = jC{v). 

Finally, we recall that if M = mi + . . . + mN is the total mass of the system, 
the center of mass of a configuration is given by the linear map G : E^ — E 

1 ^ 

i=l 

For any vector r G E wc have x • (r, . . . , r) =< MG{x)^ r >, thus it is clear that 
the space of configurations decompose as the orthogonal direct sum E'^ = A"© A, 
where X = kevG is the subspace of all centered configurations, and A the 
diagonal subspace. Wc also recall that, since the sum of all forces acting on the 
system is zero, the center of mass of any motion has constant velocity. 



1.2 The Hamilton- Jacobi equation and weak KAM solu- 
tions 

Since our mechanical system is autonomous his Hamiltonian function H is time 
independent. Therefore the natural question which arises is the existence and 
characterization of the stationary solutions of the Hamilton-Jacobi equation 
dtS + H{x,dxS) — 0, that is to say, solutions of the form S{t,x) = u{x) — ct. 
In this paper we will study the global stationary solutions, or in other words, 
functions u : E^ R satisfying the Hamilton-Jacobi equation H{x, dxu) = c, 
more explicitly the equation 



1 ^ 1 

-T — 

2 ^ mi 



du 



minii 



(1) 
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Considering the mutual distances between the bodies can be arbitrarily large, 
it is clear that equation ([T]) only has global solutions for values of the constant 
c > 0. 

A classical global solution, or global smooth solution, is a continuous function 
u : — > R whose restriction to the open and dense set of configurations 
without collisions Vt is of class , and such that ([Ij is verified at each point of 
this set. Until now, the only known smooth global solutions are solutions for 
the case N = 2, that is to say, for the Kepler problem in the space E. A direct 

1 /2 

computation shows that functions M(ri, = ± || ri — r2 || are global smooth 
solutions. 

Recently, the author has shown in P the existence of weak KAM solutions 
for a family of iV-body problems with homogeneous potentials, including the 
Newtonian case. They are global viscosity solutions of the critical equation 
H{x,dxu) = 0. The existence was obtained by application of a fixed point 
theorem to the action of the Lax-Oleinik semigroup on the convex set of viscosity 
subsolutions. 

Moreover, if F is a given symmetry group of the system, then the linear space 
of all invariant functions has nontrivial intersection with the set of subsolutions 
since clearly constant functions are invariant subsolutions. Considering that the 
semigroup preserves the set of all invariant subsolutions, which is also convex, 
the existence of invariant weak KAM solutions is also obtained (theorem 3 in 
[6]). It is evident that the application of a given isometry of the Euclidian space 
E to each component of a configuration defines a symmetry of the system. Thus 
the translation group of E, as well as the orthogonal group 0{E), acts naturally 
on the space of configurations by symmetries. Depending on the values of the 
masses, other symmetry groups can act on the space of configurations. For 
example, in the non-generic case in which the masses of some bodies are equal, 
every permutation of the position vectors of these bodies are symmetries. 

We shall remark that for a smooth, convex and superlinear Hamiltonian 
function H : T*M — M on a compact manifold M, weak KAM solutions are 
always invariant under the action of the identity component of a compact sym- 
metry group (see [5]). 

Even if the the orthogonal group is compact, there are always non-invariant 
solutions for the action of his identity component SO{E). The critical Buse- 
mann functions explicitly computed in |5] for the Kepler problem in the plane 
with equal masses, are examples of such non-invariant weak KAM solutions. 
Nevertheless, weak KAM solutions are always invariant under translations. 

Theorem 1 Every weak KAM solution of the Newtonian N-body problem in 
an Euclidean space E of dimension at least two is translation invariant. More 
precisely, given a weak KAM solution u : E^ — >■ M., we have that 

u{ri,...,rN) = u{ri +r,...,rN+r) 

for every configuration x — (ri, . . . , r^) 6 E^ and every r E E. 

The hypothesis dim£' > 2 will be only used in order to apply Marchal's the- 
orem (whose proof in the case dim£' > 3 is due to D. Ferrario and S. Terracini) 
to the calibrated curves of a given solution. This said, there is no evidence that 
suggest the theorem is not valid in the one dimensional case. 



3 



We will show that the invariance under translations of a weak KAM solution 
is equivalent to the fact that all his calibrated curves are motions with fixed 
center of mass (proposition 2]). On the other hand, a very strong property of the 
calibrated curves of a weak KAM solution is the free time minimization property, 
which we recall below. Therefore, the proof of theorem [1] will follow from lemma 
[31 where Marchal's theorem is applied. The dynamics when t — > +oo of the free 
time minimizers of the Newtonian A^-body problem is described in [T] , where is 
proved that such kind of motions are completely parabolic and asymptotic to 
some special central configurations. 

The following corollary will be easily deduced from the theorem. It shows 
as a counterpart, that the super-critical equations ([ij, that is for the values of 
c > 0, can be non-invariant under translations. 

Corollary 2 If uq : — > M is a weak KAM solution of the N-body problem, 
and r E E is an arbitrary vector, then the function 



The following section is devoted to prove the above results. It is interesting 
to remark that all the statements and proofs remains valid if we replace the 
Newtonian potential by homogeneous potential of degree a G (—2, 0). 



We start recalling some basic notions of the calculus of variations that will be 
used. The Lagrangian action of an absolutely continuous curve 7 : [a,b] E^ 



is a well defined function of 7 taking values in (0, -|-cx)]. For two given configu- 
rations x,y E E^ , the set of all absolutely continuous curves starting at x and 
arriving to y in time T > will be denoted C(x, y, T). It is well known that the 
Lagrangian action is lower semicontinuous, which implies that the infimum of 
its restriction to any set C{x,y,T) is always reached. In that follows (f){x,y,T) 
will denote the minimum of the Lagrangian action restricted to C{x,y,T). 

Let us recall Marchal's theorem on the fixed endpoint problem (cf. 
It says that minimizers of the Lagrangian action in a set C{x,y,T) not suffers 
collisions in any interior time. In other words, if an absolutely continuous curve 
7 : [a, b] — ^ E^ is such that ^(7) < A{a) for every a e C(7(a), "f{b),b - a), then 
the curve 7 avoid collisions except perhaps in the prescribed endpoints, that is 
to say, 7(t) € fl for all t € (a, b). In particular, the restriction of a such curve to 
the open interval (a, b) satisfies the equation of motion, hence it must be of class 
C°° . Furthermore, if one of its endpoints, say 7(a), is a configuration without 
collisions, then the solution can be extended to an interval of the form (a — e, b) 
for some e > 0, which implies that 7(a) exists and coincide with \imi_).^+ jit). 



u{x) = uo{x)+ < G{x),r > 



is a global viscosity solution of the Hamilton- Jacobi equation 



H{x,d^u) = 2 " ' 



2 Proofs 
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2.1 Action potential, free time minimizers and calibrated 
curves 

We will recall some results from [B] as well as some well known features of weak 
KAM theory. 

2.1.1 The action potential 

The set of all absolutely continuous curves starting at x and ending at y without 
restriction of time will be denoted C{x, y). The action potential, also called the 
Mafie critical potential in Aubry-Mather theory, is the function (p : x — 

M+ 

(j){x,y) = inf { A(7) | 7 e C(a;,y) } = |nf (/)(a;, T). 
For convenience, we will use the maximum norm in E^ ^ 
II (r-i,...,rjv) II = . max II ||^. 

i—l,...,N 

The action potential is a distance function in the space of configurations. More- 
over, it is Holder continuous in each variable with respect to any norm in E^ . 
In fact, there is a constant 77 > such that 

<P{x,y) <-n\\x-y f 

for every pair of configurations. The constant 77 only depends on the number of 
bodies N and the total mass M . Since several of the here considered notions 
are easily characterized in terms of the action potential, it will play a central 
role. 

Another important property of the action potential that we will also use, 
is that it is locally Lipschitz in 17 x 57. More precisely, if C C 51 is a compact 
set, then there is a constant A' > such that (/)(a;, y) < K \\x — y\\ whenever 
{x,y) X C. Of course the Lipschitz constant grows as well as the compact 
set C approaches some given collision configuration. 

2.1.2 Free time minimizers 

We call a free time minimizer, or semistatic curve, an absolutely continuous 
curve such that the action of each compact segment equals the action potential 
between the endpoints of the segment. In other words, a curve 7 : / E^ is a 
free time minimizer if we have ^(7 \[a^b]) — 4'{i{o)t1{^)) whenever [a, fe] C / is 
a compact subinterval. When dim£' > 2, Marchal's theorem imlpies that a free 
time minimizer defined on an interval [a, b] avoid collisions for t G (a, h). 

2.1.3 Viscosity subsolutions as dominated functions 

On the other hand we will be concerned with the notion of weak solution of ([T]) 
introduced by M. Crandall and P.-L. Lions in [3], namely the notion of viscos- 
ity solution. The classical way to define them is in terms of test functions, as 
functions which are at the same time viscosity subsolutions and viscosity super- 
solutions. It is well known that the set of viscosity subsolutions of the critical 
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Hamilton- Jacobi equation H{x, d^u) = 0, which we call the set of dominated 
functions, is exactly the convex set 

% = {u:E^ \ u{y) - u{x) < (j){x, y) for all x,y £ E'^ } . 

By the previous considerations, we know in particular that functions in H are 
Holder continuous and locally Lipschitz in n. Since £7 is a set of total measure, 
by Rademacher's theorem we also know that they are differentiable almost ev- 
erywhere. 

2.1.4 Calibrated curves 

For a given dominated function u £ %, we say that an absolutely continuous 
curve J : I ^ is calibrated for m, if for each compact subinterval [a, &] C / 
we have that u{j{b)) ~ u{j{a)) = \[a,b])- Note that a calibrated curve is 
always a free time minimizer. On the other hand, if m G H and 7 is a calibrated 
curve for u, defined on / = {a,b) and of class C^, then u is differentiable at 
7(i) for all t 6 (a, 6), and d^[t)U — C{'^{t)). It follows, using the Young-Fenchel 
inequality, that H{'~^(t),d^(t^u) = for all t e (a, 6). 

2.1.5 Characterization of weak KAM solutions 

The action of the Lax-Olcinik semigroup {Tt)t>o is well defined in [0, +00) x T-L 

by 

Ttu{x) - inf { ti{y) + 0(y, x, i) | y e } 

for t > and Tq — idu- A weak KAM solution is a fixed point of the Lax-Oleinik 
semigroup, that is to say, a function u E Ti such that Ttu = u for all t > 0. 
It is proved in |B] that such fixed points exist and that in fact, they are global 
viscosity solutions of the critical Hamilton- Jacobi equation. Moreover they are 
characterized between dominated functions as follows: u S "H is a weak KAM 
solution if and only if for each configuration x G E-^ there is a calibrated curve 
7a; : [0, +00) — !> E^ such that 7a; (0) = x. By all the above considerations, we 
deduce that if m is a weak KAM solution differentiable at a configuration a; G £7, 
then there is only one calibrated curve, since it must coincide with the motion 
generated by the initial conditions 7(0) — x and 7(0) — C^^ld^u). 

2.2 Proof of theorem [1] 

As we have said, the proof is a direct consequence of the following lemma and 
proposition. 

Lemma 3 //dim£' > 2 and x : [0, -l-oo) — >■ E^ is a free time minimizer then 
the center of mass G{x(t)) is constant. 

Proof. First, we must observe that Marchal's theorem implies x{t) G fl for 
all t > 0. Therefore x satisfies the equation of motion on (0, -f c>o), from which 
we deduce that G{x{t) has constant velocity, even at t = 0. With the notation 
Xq = x{0) and Go = G{xo), we have that G{x{t)) is of the form Go + tv for some 
vector V E E. We will prove that v — Q. 

Let y : [0, -l-oo) — t- E-^ be the internal motion associated with x, that is to 
say, the curve defined by x{t) = y{t)+6{Go+tv), where 5 : E ^ E^ is the map 
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S{r) = {r, . . . ,r). It is clear that for alH > wc have that U {x{t)) = U {y{t)) and 
G{y{t)) = 0. On the other hand, a computation shows that x-x = y-y+M \\ v ||^ 
for alH > 0. Hence we deduce that for all T > we have 

A{x\^o,T])=A{y\io,T])+lTM\\vf. 

Since a; is a free time minimizer, and cp verify the triangle inequality, we also 
obtain 

A {x |[o,T]) = (t>{xo, xiT)) < (j) {xq, xq + T5{v)) +(l){xo+ TS{v), x{T)) . 

Let yT be the curve obtained translating the restriction of y to the interval [0, T] 
in such a way that the center of mass of yT becomes fixed at G{x(T)). Thus, 
yx is the curve defined for t G [0, T] by yrit) = y{t) + S{Go + Tv). Both curves 
have the same action, but now we have 7t G C{xo + TS{v), x{T)). Therefore, 

0(x(O)+T<5(«),x(T)) < A(j/|[o,T]). 
We use now the Holder estimation for the action potential, and we get the bound 

(f> {x{0), x{0) + T5{v)) < r] \\ T5{v) ||^/^ = T^/^ r] || v ||^/^ . 
We can conclude that the inequality 

lTM\\vf < ri/2^||^||V2 

holds for all T > 0, but this is possible only in case v = 0. □ 

Proposition 4 A weak KAM solution u : — >■ M is invariant under transla- 
tions if and only if all calibrated curves have constant center of mass. 

Proof. We start showing the necessity of the condition. Let u be a weak 
KAM solution which is invariant under translations. At each point of differ- 
entiability X G we must have have A c kerrf^w- Let x G E^ be any 
configuration, and let '■ [0, +00) be a calibrated curve for u starting at x. We 
know that u is differentiable at "fx{i) for alH > and that d^^(^t) = ^{ixit)), 
that is to say, dxu{w) = w-jxit) for all w G E^ . It follows that '^x{t) G A-"- = X 
for t > 0, hence that G{jx{t)) is constant. 

Let now u be a weak KAM solution such that all calibrated curves have 
constant center of mass. Let X G £^ be some given configuration, and choose 
a non zero vector r G E. We will prove that u{x) = u{x + S{r)), where 6{r) = 
(r, . . . , r) G A. 

Since u is continuous and the set Cl is dense and translation invariant, it 
suffices to consider configurations x £ Cl. Let 

B{e,r) = {zeE'' \z-d{r) = 0, \\z\\<e} 

be the closed ball of radius e > in the orthogonal complement of S{r) G A. 
Since a; G O, we can assume that e is so small that the set 

B = {x + z\ze B{e,r)} 
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is contained in n, which implies that also the compact cylinder 

C={y + tS{r)\ijeB,te [0,1]} 

is contained in il. Let us call C" C C the set of points where u is difFerentiable. 
Recall that w is a dominated function, hence it is a Lipschitz function on C, 
and C has total measure in C. For q g C", let -fg : [0, +00) be the only 

calibrated curve for u with 7^(0) = q. We know that dqu{w) = w ■ jq{0) for all 
■w S . By using the hypothesis we deduce that 7g(0) £ X = kerG = A-*-, 
and therefore that A C kerd^u for all q € C . 

We consider now the Lipschitz continuous function / : S x [0, 1] — ?> K defined 
by f{y,t) — u{y + t5{r)) — u{y). It is clear that dtf is zero almost everywhere. 
Therefore, applying Fubini's theorem dtf on a cylider of the form A x [0, 1], 
where A C i? is an open set, we obtain 

/ /(y,l)dy = 0. 

J A 

Since the open subset A d B \s arbitrary, we must have f{y,l) = for all 
y & B. Hence, we have /(x, 1) = 0, that is to say, u{x + S{r)) ~ u{x). □ 

2.3 Viscosity solutions of the supercritical equations 

Proof of corollary\2\ Let uq be a weak KAM solution and r £ E. We will use 
the fact that uq is a global solution of H{x,dxu) = in the viscosity sense. As 
before, we call S : E ^ A the trivial isomorphism. Thus function Ur can be 
written as Ur = mq + Or, where is the linear function ar{x) = x ■ S(r). 

Let Xo be a given configuration, and suppose that (p G C^{E^) is such 
that if — Ur has a local minimum at xq. Since uq is a viscosity solution of 
H{x,dxu) — and we have that (ip — a^) — uq has a local minimum at xq, we 
have H{xo, d^giv ~ O'r)) < 0. 

On the other hand, by theorem [1] we know that uq is translation invariant, 
hence A C ]ieTdxo{^ — a-r). Since dx^ar = Or, it is clear that X C kerdxaO-r and 
that I dx^ar \ = M* \\r\\ . Therefore we deduce that 

I dxa^P f = I dxo{(p - Qr) f + M* \\ r ||^ 

and we conclude that 

H{xo,dxo^) = H{xo, dx.iifi - ar)) + ^M* \\ r f < ^M* \\r\\\ 

proving that Ur is a viscosity subsolution of the supercritical equation at a;o. 
With a symmetric argument it can be proved that Ur is a viscosity supersolution 
of the same equation. □ 
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